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Origins of Taylor's power law for fluctuation scaling in complex systems 
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Taylor's fluctuation scaling (FS) has been observed in many natural and man-made systems 
revealing an amazing universality of the law. Here we give a reliable explanation for the origins and 
abundance of Taylor's FS in different complex systems. The universality of our approach is validated 
against real world data ranging from bird and insect populations through human chromosomes and 
traffic intensity in transportation networks to stock market dynamics. Using fundamental principles 
of statistical physics (both equilibrium and non-equilibrium) we prove that Taylor's law resuhs from 
the well-defined number of states (NoS) of a system characterized by the same value of a macroscopic 
parameter (i.e., the number of birds observed in a given area, traffic intensity measured as a number 
of cars passing trough a given observation point or daily activity in the stock market measured in 
millions of dollars). 

PACS numbers: 89.75.-k, 89. 75. Da, 05.40.-a 



I. INTRODUCTION AND MOTIVATION 

In ecology, Taylor's power law [l| (or the law of the 
mean) states that the mean, {N) , and the variance, cr^ = 
(iV^) — (iV)^, characterizing the number of population 
representatives are related by power law. 



(7V2) - (7V)2 = a{N)\ 



(1) 



with the characteristic exponent, 6, describing effects of 
heterogeneity in spatial or temporal patterns of the fre- 
quency distribution. The value of h is usually in the range 
of 1 to 3. For comparison with the Poisson distribution, 
b — \, the parameter h > 1 corresponds to clustering 
(aggregation), whereas 6 < 1 may be interpreted as or- 
dering. 

As we have already stated, the FS described by Eq. ([T]) 
has been noted in a variety of natural and man-made sys- 
tems, and the universality of Taylor's law is now widely 
recognized 0-lE| • To emphasize that the generality of our 
approach to this law is commensurate with the universal- 
ity of the law itself, we do not concentrate on a specific 
system in this paper. Later in the text, the macroscopic 
quantity N simultaneously stands for the number of birds 
observed in a given area, the number of cavities and corn 
borers found on one plant, the number of gene struc- 
tures located in equal-sized non-overlapping bins that 
span the whole chromosome, the daily traded value of 
a given stock at the New York Stock Exchange (NYSE) , 
and the number of cars passing through a certain obser- 
vation point in a given time period. 

When the data over which the averages ([T]) are taken 
have a temporal structure (like in the case of traffic) , the 
Taylor's law is called temporal FS, otherwise (like in the 
case of chromosomes) we term it ensemble FS [Ij. We 
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show, that given our approach one does not need to in- 
voke any stochastic models to explain phenomena such 
as aggregation effects in different populations, and traffic 
jams. Regardless of character of the scaling (ensemble or 
temporal) the collective phenomena manifesting them- 
selves by Taylor's law simply result from properties of 
the phase space underlying the considered systems. 

The only assumption made concerning the systems 
obeying the Taylor's law is that they are either in equi- 
librium or in non-equilibrium steady states. The equi- 
librium formalism described in the paper is suited for 
systems with ensemble FS that are characterized by sta- 
ble, time-independent probability distributions. Systems 
with temporal FS described by time series with a sta- 
ble (i.e. constant in time) mean value and variance, are 
treated as non-equilibrium systems in steady states. In 
the paper, the case of ensemble FS is represented by the 
data on spatial distribution of larval populations of the 
European corn borer and the data on spatial distribu- 
tion of gene structures on human chromosome 7, while 
the case of temporal FS is represented by traffic in Min- 
nesota and the steady-state dynamics of North American 
avifauna and NYSE. 



II. DATA DESCRIPTION 

Here, we briefly describe the datasets used to validate 
our theoretical approach. We also draw readers attention 
to some known problems (or new observations), which are 
somehow related to Taylors law and can be explained in 
terms of the theoretical framework. 

European corn borer. A few equilibrium, frequency 
distributions describing larval populations of this pest 
(Pyrausta nubilalis) have been published in a 1957 paper 
in Biometrics '6*1 . The paper has brought up an intrigu- 
ing (and until now unsolved) issue of what kind of fre- 
quency distribution should be used to describe different 
populations. In the study, three areas, each consisted of 
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FIG. 1: Taylor's fluctuation scaling for systems analyzed in the paper, a, European corn borers {Pyrausta nubilalis). 
b, Gene structures in human chromosome 7. c, Daily turnovers of stocks traded on the New York Stock Exchange (NYSE), d, 
Mallards {Anas platyrhynchos) and Blue Jays ( Cyanocitta cristata) in all routes populated by these species in 1966 — 2007. e, 
Savannah Sparrows {Passerculus sandwichensis) in all routes populated by this species in 1966 — 2007. f, Savannah Sparrows 
counted in routes corresponding to two difFcrcnt physiographic conditions located in the Drift Prairie and the Southern Rockies, 
g, TrafBc intensity as measured by all automatic traffic recorders in Minnesota in 2007. h, Daily fluctuations in traffic measured 
by a single recorder (ATR no. 222) in 2002 — 2007. In the figure, fiiU points correspond to raw data (i.e., mean and variance 
calculated according to the description given in the text and in the Table I), the open symbols express logarithmic binning of 
the data, and the solid lines represent their linear fits (in log- log scale). In the case of birds, the key observation is that the 
characteristic parameter 6 may differ not only among the species but also within the same species if one takes into account 
the physiographic stratification characterizing living conditions in different areas (for that reason for further analysis we have 
selected Blue Jays as a species that weakly depends on the physiographic stratification). A similar comment is true for the traffic 
intensity. In the latter case, one can see that although the whole data obey Taylor's fiuctuation scaling with the parameter 
h = 1.43, in reality the data are very heterogeneous. The parameter b characterizing traffic intensity recorded by a single ATR 
may change dramatically from hour to hour. 
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Panel in Fig. 1 


No. of 
points 

panel 


No. of data from 
which each point in 
the panel (i.e. mean 

and variance) has 
been calculated 


Individual data 
represents 


Comments 


la 


3 


1296 


Total number of corn 
borers living on a given 
plant 


Each point in the panel represents mean and variance for one 
of three areas each consisted of approximately 3 acre, chosen 
from the same corn field. 


lb 


8 


40 - 3168 (variable) 


Total number of genes 
positions of which start 
in a given bin 


C'hi omosome T sp ans about 159 million H) N' A. building bio cks 
(base pairs) and contains more that 2100 genes (including 
pseudogenes). The chromosome has been divided into equal- 
sized non-overlapping bins that spanned the physical length of 
the chromosome. Number of data from which mean and vari- 
ance have been calculated depends on the width of the bins 
(from 5 X 10^ to 4 X 10^). 


Ic 


4728 


10 


Daily turnover of a 
given stock 


Daily turnovers of 2364 stocks were retrieved over two periods 
of ten consecutive trading days. 


Id (Mallard) 


3310 


5 - 210 (variable) 


Total number of indi- 
viduals of the species 
recorded in a given part 
of a given route in a 
given year 


Each route has been divided into 5 parts (stops 1-10, 11-20, 
21-30, 31-40, and 41-50). Maximal value in column 3 follows 
from the calculation: 210 — (5 stops on a route) X (41 years). 
Because some species were not recorded each year this quantity 
can be lower. The value in column 2 is just a number of routes. 


Id (Blue Jay) 


3017 


le 


2428 


If (Prairie) 


73 


ii ( IxocKies ) 


51 


Ig 


41472 


20 - 23 (variable) 


Total number of cars 
passing a given ATR in 
a given hour in one di- 
rection in a given week- 
day of a month of a 
given year 


The value in column 3 is equal to the number of weekdays in 
a month (so it can vary from 20 days to 23 days depending on 
the month). The value in column 2 follows from the calcula- 
tion: 41472 ^ (72 ATRs) x (2 directions) x (24 hours) x (12 
months) X (1 year) 


Ih (lain-7am) 


504 


Value in column 2 follows from the calculation: 504 — (1 A'l'H) 
X (l direction) X (7 hours) x (12 months ) x (6 years) 


Ih (8am-6pm) 


792 


Value in column 2 follows from the calculation: 792 — (1 ATR) 
X (1 direction) X (11 hours) X (12 months) X (6 years) 


Ih (7pm-12am) 


432 


Value in column 2 follows from the calculation: 432 — (1 ATR) 
X (1 direction) X (6 hours) x (12 months) x (6 years) 



TABLE I: Detailed description of data used to compile Figure 1. 



approximately 3 acre, chosen from the same corn field, 
have been investigated. For each area a distribution of 
larvae per plant has been obtained. In order to fit the 
gathered data its authors considered three different com- 
pound Poisson distributions: the negative binomial, the 
Neyman type A, and a distribution they called Poisson 
binomial. Before and a long time after its publication, 
the distributions described in the paper have been used to 
characterize different diversity patterns in population bi- 
ology [S]. Later in this article, we derive a correct formula 
for the frequency distribution describing the populations 
that reveal Taylors law of the mean. 

Human genome. Data on physical distribution of gene 
structures on human chromosome 7 were retrieved from 
a database provided by the Chromosome 7 Annotation 
Project [13] • Previous analysis of the data has demon- 
strated that the density of gene structures within the 
chromosome is heterogeneous . It has been also shown 
that the number of such structures contained within a se- 
quence of equal-sized non-overlapping bins that span the 
physical length of the chromosome fulfills Taylors fiuctu- 
ation scaling with the characteristic exponent, b — 1.61 
0. This observation has been recognized as a quantita- 
tive test confirming the presence of gene clustering within 
the human genome. In this article, we show that Taylors 
law with 6 > 1 is always due to collective phenomena 
such as clustering effects. The case of 5 = 1 corresponds 



to complete randomness described by the Poisson distri- 
bution. 

Traffic in transportation networks. Hourly numbers of 
cars passing through observation points located on inter- 
states, trunk highways, county state-aid highways, and 
municipal state-aid streets at various locations through- 
out Minnesota were retrieved from the Minnesota De- 
partment of Transportation [45]. The traffic intensity 
had been recorded by 72 automatic traffic recorders 
(ATR) from 2002 to 2007. The datasets, from which 
mean and variance were calculated, include the number 
of cars observed by a single recorder and passing in only 
one direction at a given hour in all weekdays of a month 
Q. We have checked that such time series have a sta- 
tionary character (i.e. traffic is homogeneous within the 
considered periods). In Figure [TJ one can see that al- 
though the whole data obey Taylors fluctuation scaling 
with the parameter b = 1.43, in reality, the data are very 
heterogeneous and the parameter, 6, characterizing the 
data recorded by a single ATR may change dramatically 
from hour to hour. 

North American avifauna. The data were retrieved 
from the North American Breeding Bird Survey [i^ . The 
survey has collected annual abundances over a 40-year 
period for over 600 bird species in more than 3000 ob- 
servation routes. In our analysis, we have concentrated 
on single species detected in consecutive years on sepa- 
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rate routes. Given such a time series (i.e. the number 
of birds of a given species counted in a given route in 
consecutive years), one can calculate its mean and vari- 
ance (it has been proved elsewhere that the considered 
time series have a stable temporal structure, see We 
have done so for three rather abundant species represent- 
ing different bird families: Mallard (Anas platyrhynchos, 
Family Anatidae), Blue Jay {Cyanocitta cristata, Fam- 
ily Corvidae), and Savannah Sparrow (Passerculus sand- 
wichensis, Family Emberizidae). Figure [T] presents the 
empirical Taylors law for these species. The clue obser- 
vation is that the characteristic parameter, b (see Eq.[T|), 
may differ not only among the species but also within the 
same species if one takes into account the physiographic 
stratification characterizing living conditions at different 
routes (see for speculations about the meaning and 
the value of the parameter b). 

New York Stock Exchange (NYSE). Daily turnovers of 
2364 stocks were retrieved over two periods of ten consec- 
utive trading days (May 8 - May 21, 2008 and May 22 - 
June 5, 2008) from Yahoo Finance Stock Research Center 
(47t . The periods were selected due to relative stability 
of the NYSE Composite Index and the trading activity 
in the market. Mean and variance were separately cal- 
culated for each period and each stock. Taylors scaling 
exponent of the system, b = 1.87, was quite large. We 
show that the large value of the parameter, b, and Levy 
flights, which are commonly observed in stock market dy- 
namics [H, [l^, are related to each other, as they result 
from non-trivial properties of the phase space underlying 
the considered systems. 

Detailed description of how the data were processed is 
given in Table I and in Appendix [X] 



III. ENSEMBLE FLUCTUATION SCALING 

To explain the origins of Taylor's power law for ensem- 
ble FS in different complex systems, we start with equi- 
librium statistical physics. It is well known that the dis- 
tribution 7'(f2;/i) constrained to yield the average value 
of the parameter N is given by 14] 



(2) 



where n stands for the external field coupled to N that 
imposes a given value of (N) and F{^) represents the so- 
called free energy of the considered system that encodes 
properties of the system in equilibrium. 

Here, the crucial point to understand is that the Greek 
letter ft in Eq. ([2]), refers to the so-called microstate 
of the considered system, whereas researchers studying 
real-world systems are usually interested in macrostates 
and the corresponding macroscopic quantities. Accord- 
ingly, in this investigation, instead of the distribution 
V{fl;fi), we seek the expression for the frequency dis- 
tribution P{N; /i) (sometimes we write P{N) instead of 
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FIG. 2: Comparison of real-world and theoretical fre- 
quency distributions P{N) characterizing the consid- 
ered systems for the given values of (A'^) (left column), 
and quotients of two distributions fitted by the expo- 
nential function predicted by Eq. (|18p (right column). 
The figure consists of 8 panels arranged in 4 rows. Each row 
corresponds to a diflterent dataset and provides information 
on (for detailed information on the data used see Section |ll] 
and Appendix [X] detailed description of how fitting of ex- 
perimental frequency distributions has been done is given in 
Appendix IB 4|l ; a, b, European corn borer frequency pat- 
terns; c, d, Intensity of car traffic in Minnesota, as measured 
by a single recorder (ATR no. 222) during night hours (1 
a.m. - 7 a.m.) in 2002 — 2007; e, f, Blue Jay (Cyanocitta 
cristata) abundance in North America in 1966 — 2007; h, g, 
Daily activity at the New York Stock Exchange (NYSE) in 
the period May 8 - June 5, 2008. Some discrepancies visible 
in the last case may be due to the problem with selection of 
more homogeneous subset of analyzed stocks. 



P{N; /i)) that characterizes macroscopic states of differ- 
ent systems fulfilling Taylor's law. The essence of our 
approach lies in a rather trivial observation that the two 
introduced distributions are related to each other by a 
simple expression. 



P{N■,^l)^g{N)V{n;^l), 



(3) 



5 



where g{N) is the announced NoS, which gives the num- 
ber of microstates 51 having the same value of the macro- 
scopic parameter N . In the foUowing, we show that in 
the systems fulfiUing Taylor's FS, the function g{N) has 
a well-defined form, and this immediately allows us to in- 
fer the origin of Taylor's law to understand the meaning 
of Taylor's characteristic parameter, b. 

Provided that N is nonnegative and discrete, the nor- 
malizing factor in Eq. ([2|), e^^^\ is just the Z-transform 
ofNoS,g(iV), 



(4) 



N=0 



Using the exponential formula known from combinatorial 
mathematics [i5[, the left-hand side of Eq. Q can be 
written as 



^=1 



g/o 

-^-B7v(/i, /2, • ■ • , In) 



f^'\ (5) 



where i?Ar(/i, /jv) is the A^th complete Bell polyno- 
mial whereas /„ represents the coefficient of the nth term 
in the MacLaurin expansion of the free energy. 



OO p 



(6) 



ri=0 



In accordance with Eq. ([5]) , the right-hand side of Eq. ^ 
can be written in the following form 



E 



E(-i)" 



JV=0 



ji=0 



.N=0 



(7) 



Although the expressions in ([5|) and (O seem to be rather 
complicated, their theoretical and real-world interpreta- 
tion is very simple. 

Let us concentrate on mathematical issues. First, the 
formula within the square brackets in Eq. ^ is the Pois- 
son transform, 



Gin) = P.T. [GiN),n] = ^ G(7V)e-^7V"/n!, (8) 
of the function 



N=0 



G{N) = e^g{N) 



(9) 



When working with the Poisson transform, it is impor- 
tant to understand how the transform acts on an arbi- 
trary function. Simplifying, one could say that the trans- 
form adds Poissonian fluctuations to the function. For 
that reason, Q (n) looks like a fuzzy image of the original 
G{N), and it is often reasonable to assume that (for more 
details see Appendix IB ip 



(10) 



g(n) =P.T. [G(iV),n] ~G(n). 



The last property of the transform turns out to be very 
useful in the inverse problems, in which one has to calcu- 
late the original function provided that its Poisson trans- 
form is known. Such a problem arises in our derivations. 
Comparing Eqs. ([5]) and ([7]) for n > 1 one gets 



5{n) 



(-l)"li?„(/i,/2,...,/„). (11) 
n] 



The closed form of the inverse Poisson transform for Q (n) 
and the closed form of g{N) only exist in a few cases. In 
the general case, for iV ;» 1, one can use the following 
approximation 



-N 



g{N)^e-"giN)<x—Bj,ihj2,...jN). (12) 

Eq. (IT^ is the main theoretical result of this article. It 
has a lot of in common with the famous Mayer's diagram- 
matic expansions for imperfect gas (l^ . [la | . Further in 
this paper we show that the derived NoS function has a 
very intuitive form. 

To apply the derived formula to systems with ensemble 
FS one has to know all the parameters /i, /2, . . . , /at, i.e. 
one has to find F{fi) describing systems obeying Taylor's 
law. In order to do it we exploit fluctuation-dissipation 
relation 



d{N) _ d'^Fifi) 



(13) 



which states that fluctuations of the parameter N are 
proportional to susceptibility of the parameter to its con- 
jugate field /i. Comparing right-hand sides of Eqs. (fT3|) 
and (HI one obtains differential equation for (A''), i.e., 
—d{N)/dfi = a{N)''. Solving this equation, with the 
reasonable assumption of nonnegative variance, one gets 



(N) 



for 5=1 



((6- l)a/i-|- X)i/(i-'') for 6>1. 



(14) 



Next, having and again exploiting (|13p . i.e., solving 
dF/dfi = —{N), one also finds the formula for the free 
energy, F{fi), 

( 1 (N)'^'^^''') +Y for b>l 

X and Y represent integration constants in Eqs. 
and dTSl) . Coefficients in the MacLaurin expansion ([5]) 
of the free energy F(fi) are given by (for derivation see 
Appendix |B2)) 



(16) 



where 



c„ = (-l)"a"-i l[[{n - 2)6 - (n - 3)]. (17) 



j=2 



Detailed calculations for the Poisson distribution 
P{N; fi) ([3]) are shown in Appendix IB 31 
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FIG. 3; Experimental verification of Eq. (flOl). Differ- 
ent symbols placed in the graph represent different real- world 
systems that show Taylor's fluctuation scaling. Note that the 
symbols cover the solid theoretical line to an impressive extent 
of nine orders of magnitude of Afi. 



Now, before delving into an interpretation of the de- 
rived expressions, it would be valuable to convince the 
reader that our approach does really account for the be- 
havior of real- world systems with ensemble FS. To do so, 
we performed three quantitative tests on experimental 
data. The tests clearly show that Taylor's law results 
from the unique NoS underlying the considered systems. 

The first test consists of a direct comparison between 
experimental frequency distributions and theoretical dis- 
tributions given by Eq. ([3]). The second test shows that 
the quotient of two frequency distributions correspond- 
ing to different average values of (N) is an exponential 
function of A^, 



Jf^2-fJ.i)N 



(18) 



The third test follows from Eq. and indirectly refers 
toEq. ((T8)) . Namely, transforming the expression for (N), 
one can show that the experimental data should satisfy 
the identity 



- Ml = 



{N2 



a(l - b) 



= h{{N,),{N,)). (19) 



Figures [2] and |3] present results of the three tests applied 
to systems with ensemble FS: European corn borers and 
human genome. Note the excellent agreement between 
the data and our theoretical predictions. 



IV. TEMPORAL FLUCTUATION SCALING 

Now, the question is: Does the explanation of ori- 
pn of ensemble FS in equilibrium systems may help to 



understand temporal FS characterizing systems in non- 
equilibrium steady states? The answer is affirmative. Re- 
cently, it was shown that phase space probability distri- 
bution describing such systems has an exponential form 
[l3[il, cf. Eq. ©, 



(20) 



where Ni is the so-called first moment of an additive vari- 
able N that the considered system may exchange with the 
reservoir. The parameters fi and fj,i stand for external 
fields coupled to N and A^i respectively, whereas M^™"' 
is called the mirror work. Finally, F(ll, ni) represents 
nonequilibrium free energy (see also 19]). 

The formula (|^n|) derived by Phil Attard generalizes 
the Boltzmann distribution, Eq. to nonequilibrium 
systems. To account for temporal FS described by 
Eq. dll), we use the formula to calculate frequency dis- 
tribution P{N] jjL, Hi). Similarly as in the case of equi- 
librium statistical physics, to get P{N; h, fii) one has to 
sum Eq. (|20p over all microstates, that fulfill the con- 
dition N{Cl*) = const, i.e. 



P(7V;/i,Aii) = ^7'(r!;/i,Aii) = 



-fiN 



(21) 

where g*{N;fii) = J^n* g-t^iiNi+wrn may be termed 
as the weighted NoS. Finally, averaging Eq. (|211) over 
different values of the parameter ni one can reduce it to 
the formula that is equivalent to Eq. (jS]) 



{PiN■,^i,^l^)) =PiN;^l,^,l) (22) 



PiN;^i) 



where the last transformation in (j22p is valid under the 
generalized mean value theorem and fi* is a certain value 
of the parameter /ii. 

The reasoning behind the algebra is that the formalism 
applied to ensemble FS may be adopted to explain tem- 
poral fluctuations. Results of three tests (as described 
earlier) applied to systems with temporal FS (i.e. traf- 
fic intensity in Minnesota, dynamics of North American 
breeding bird populations and NYSE) shown in Figs. [5] 
and [3] certify the statement. 



V. NUMBER OF STATES IN REAL- WORLD 
SYSTEMS 

In the following, we concentrate on the formula for 
the number of states, g{N) p^ . Although the formula 
seems rather complicated, its theoretical and real-world 
interpretation is very simple. In general, the A'^th com- 
plete Bell polynomial, i3jv(/i, f2: - ■ ■ , Jn), describes the 
number of disjoint partitions of a set of size N into an 
arbitrary number of subsets [l^, [2l| . The parameters fi 
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with i = 1,2, . . . , N apply to subsets of size i and play an 
important role in a description of the partitions. For ex- 
ample, if all the parameters fi, f2, ■ ■ ■ , /n have the same 
value, then there is no preference for the size of sub- 
sets. The resulting partitions correspond to a random 
distribution of elements, and the frequency distribution 
describing that system, P{N) ([3|), is Poissonian (for a de- 
tailed derivation see Appendix IB 3p . On the other hand, 
in the extreme case of fi ^ fj for all j i, the Bell 
polynomial gives the number of such partitions in which 
there is a strong preference for subsets of size i. 

Given the meaning of the complete Bell polynomials, 
an interpretation of the number of states (fT2|) underlying 
real-world systems with Taylor's fluctuation scaling fol- 
lows immediately. For example, in the case of the bird 
population of size N, the number of states, g{N), is pro- 
portional to the number of different partitions of N birds 
into subpopulations of arbitrary size. The number and 
the size of subpopulations are encoded in both the free 
energy of the system, F{^), and the corresponding pa- 
rameters, fi. The analogous interpretation of the NoS 
applies to every other animal, insect, and plant popula- 
tion. 

In the case of stock market dynamics, traffic in trans- 
portation networks, and other systems driven by human 
activity, the number of states in (|12p has a similar in- 
terpretation. In order to show the analogy, let us con- 
centrate on the number of cars (e.g., N = 60) passing 
through a given observation point in a given time period 
(e.g., T = 60 minutes). First, note that the parameters 
N and T do not tell the whole story about the traffic. 
The same values of TV and T may result from homoge- 
neous traffic (e.g., on average one car per one minute) or 
inhomogeneous traffic (e.g., all cars counted within the 
first five minutes). According to our approach, all mi- 
crostates il with the same N are equiprobable as shown 
in This does not mean, however, that all the states 
are possible from the point of view of the considered sys- 
tem. For example, if the observation point is located on 
a very busy roadway, then states of the roadway with no 
cars running are unlikely. 

In fact, the problem of car traffic may be simplified 
to the problem of N balls in T boxes (boxes may repre- 
sent minutes with balls corresponding to cars). In this 
notation, the meaning of the number of states given by 
the complete Bell polynomials is clear. The number of 
states, g{N), corresponds to the number of different par- 
titions of N balls (cars) into T boxes (minutes). The 
parameters /i, /2, . . . , /jv indicate which partitions are 
reasonable and likely from the point of view of the con- 
sidered system. If the balls are noninteracting, all the 
parameters fi have the same value. On the other hand, 
when the spectrum of fi is not uniform, there must exist 
interaction between the balls (cars) resulting in different 
collective phenomena, e.g., clustering seen as a traffic 
jam. 

To complete the discussion of the number of states 
given by Eq. ([T^ . one should mention the so-called par- 



tial BeU polynomials, BN^k{fi,f2, • ■ • , fw-k+i) [Hi. The 
polynomials describe the number of partitions of a set 
of size N in which exactly k subsets are considered, 

-BAr(/l, /2, • ■ • , /w) = -Sjv,fc(/i, /2, ■ • ■ , /w-fc+i)- 

The decomposition of the complete Bell polynomials into 
partial polynomials allows one to analyse how different 
partitions contribute to the number of states for different 
values of the parameter b. 

In Figure m one can see such a decomposition for dif- 
ferent values of N and 1 < 6 < 3. It is remarkable 
that the preferred number of subsets, (k), and the pre- 
ferred size of these subsets, N/ (k), have well-defined val- 
ues, and these values strongly depend on the parameter 
b of Taylor's power law ([1]). The average number of sub- 
sets is a decreasing function of b and, similarly, the av- 
erage size of subsets is an increasing function of b. This 
means that larger values of the parameter characterise 
stronger collective phenomena manifesting either in ag- 
gregations of individuals in different populations [l|, llll ] , 
traffic jams [121, and Levy flights in stock market dynam- 
ics [ill . This supports the scientific message of this arti- 
cle; that Taylor's fluctuation scaling is due to the number 
of states underlying the considered systems. The number 
of states has a built-in susceptibility of the system to col- 
lective effects, the strength of which depends on the value 
of Taylor's parameter, b. It also means that one need not 
invoke any stochastic models to explain these phenom- 
ena. In fact, there may exists a number of stochastic 
processes deflned in the phase space with NoS given by 
Eqs. and p5|) . Of course, all the processes will result 
in Taylor's fluctuation scaling. 



VI. CONCLUDING REMARKS 

In closing, we note that our approach differs crucially 
from previous work on Taylor's power law [llj, [23l - [30| . It 
is significant that the approach does not invalidate the 
previous models. On the contrary, our approach may al- 
low identification of the models' important features that 
account for the law. We show that Taylor's fluctuation 
scaling results from the ubiquitous second law of ther- 
modynamics (here called the maximum entropy princi- 
ple) and the number of states (a concept borrowed from 
physics). We anticipate that our formalism will provide a 
quantitative basis for formulating a new theory of popu- 
lations, communities, and ecosystems [^, a theory based 
on concepts of the number of states. In epidemiology 
and medicine, our explanation of Taylor's law may be 
helpful in accounting for such observations as clustering 
in human sexual contact in HIV transmission [31], epi- 
demic outbreaks [33 | , or high variability of cancer statis- 
tics within the human population (33| . We believe that 
the approach may also help in understanding the variabil- 
ity in organ cell numbers for a variety of organisms [34| 
and the clustering of genes on human chromosomes [8[. 
In all the listed cases, Taylor's law has been recognized 
as an intrinsic feature of the considered systems, and the 
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FIG. 4: Interpretation of the number of states function g{N) given by Eq. (|12|l . Decomposition of the complete Bell 
polynomials into partial Bell polynomials for different values of Taylor's parameter b. We have already shown that different 
values of b result from different free energies F{fi) (|15[) characterizing the considered systems, which, in turn, provide different 
sets of the coefficients fi, f2, ■ ■ ■ , /n ■ Graphs a and b show that the whole set of partitions for a given b is dominated by 
partitions that consist of a well-defined number of subsets, k (here, A'^ = 50). This number decreases with b leading to an 
increase in the average size of these subsets, N/{k) (cf. graph c). The decomposition analysis clearly shows that higher values 
of b correspond to aggregation effects. The elements of the original set of size aggregate into subsets whose sizes increase 
with b. 



new interpretation of the law offered here may help clar- 
ify the underlying dynamics. 
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Appendix A: Data processing 

We would like to stress that although there are plenty 
of real data in which Taylors fluctuation scaling can be 
observed, most of the data require careful processing be- 
fore they can be used in our analysis. The two basic 
problems that we encountered in trying to validate our 
theoretical approach in real-world systems were related 
to: heterogeneity of the considered data and/or small 
amount of the data. 

The first problem is partially illustrated in Fig.[T] The 
essence of the problem is that, although the bulk of the 
considered data fulfills Taylors law with a given parame- 
ter b, in fact, the data may consist of a certain number of 
subsets i = 1, 2, . . . (e.g., representing logical subsystems 
of the considered system), each of which is characterized 
by its own characteristic parameter, bi. Careless selection 
of data may lead to misleading results. The quantitative 
tests basing on Eqs. ([3]), (fTO)) and described in the 
article are peculiarly sensitive to such a careless attitude 
to data. Therefore, in most cases, data standardization is 
required. Unfortunately, such a standardization may sig- 
nificantly decrease the amount of data, ft may even lead 



to questions on statistical significance of the preformed 
analysis. 

The problem of statistical significance is the more im- 
portant because, although Taylors law ([1]) operates on 
macroscopic parameters (i.e., on the mean, (N), and the 
variance, cr^ — {N'^) — {N)"^, of the frequency distri- 
bution, P{N)), from both experimental and theoretical 
points of view the basic observable in our approach is the 
distribution, P{N), itself. To prove validity of our ap- 
proach, one has to operate on large datasets, which pro- 
vide smooth experimental frequency distributions, P{N), 
and are also homogeneous in the sense of the parameter, 
b. 

For these reasons, the actual verification of our theo- 
retical approach shown in Figs. [5] and [3] m the main body 
of the article was, inter alia, based on the most abundant 
bird species representing the North American avifauna 
and on the frequency distributions describing larval pop- 
ulations of European corn borer published in a classical 
paper [Q] . In analyzing traffic intensity in Minnesota, we 
have concentrated on only one (from among 72 others) 
traffic recorder in night hours. To increase the amount 
of data, we have analyzed traffic over a very long time 
period of 5 years (2002 - 2007). Analysis of such a long 
time period was possible due to relative stability of traf- 
fic measured by the considered recorder during this time 
period. A similar procedure of data selection has been 
applied to the NYSE data. However, although we ex- 
pected the system to be very heterogeneous in the sense 
of Taylors parameter, 5, we did not observe any of these 
effects. Contrary to transportation network, our only 
concern with the stock market was the highly nonsta- 
tionary character of the data. Therefore, a more detailed 
analysis of NYSE was based on daily turnovers of all 2364 
stocks quoted in this market in a rather short time pe- 
riod 2x10 consecutive trading days (i.e. May 8 - May 
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o 



logio<N) 

FIG. 5: (Color online). Schematic description of the 
methodology used in real data analysis. Detailed de- 
scription is given in tlie text. 



21, 2008 and May 22 June 5, 2008). The periods corre- 
spond to relative stability on both trading activity in the 
market and the NYSE Composite Index 39. 

The above description of data selection together with 
the general description of the considered datasets should 
make possible independent reconstruction of Fig. [TJ Be- 
low, we describe our method to obtain the smooth exper- 
imental frequency distributions, P{N), shown in Fig. [51 

We were interested in smooth frequency distributions, 
P{N), possessing the given mean, (iV), and the corre- 
sponding variance, ajf = {N'^) — (N)'^. Having the mean- 
variance graphs, as shown in Fig. [TJ the smooth distri- 
butions characterizing single points in these graphs can 
be only observed in two cases of distributions describ- 
ing larval populations of European corn borer and gene 
structures in the human chromosome 7. In the remain- 
ing cases (including bird populations, NYSE, and trans- 
portation network) , distributions corresponding to single 
points on the mean-variance graph are very noisy. To 
increase the amount of data from which a distribution 
is made, we have assumed that the neighboring points in 
the mean- variance graph result from similar environmen- 
tal conditions (i.e., from similar values of the parameter 
fi). In this way, a single smooth distribution possess- 
ing the given values of (N) and a% has been prepared 
as a simple sum of all the component distributions cor- 
responding to single points in the mean-variance graph 
and meeting the following conditions: 



log 



10 



and 



log 



(N) 



'N,i 



10 ■ 



'N 



< logio d 



< logio d 



(Al) 



(A2) 



Interpretation of the parameter d is easy. It describes 
a linear size of the square in the log-log plot in the 



mean- variance graph with the central point of the square, 
[{N),aj^], placed in the solid line corresponding to the 
empirical Taylors law (cf. Fig. [5]). If d is chosen to be 
too small, i.e., d — > 1, then the square shrinks to the 
single point, [(iV),cr^]. It leads to poor data samphng 
and noisy distribution, P(N). On the other hand, when 
the parameter d is taken too large, the resulting distri- 
bution, P(iV), is made of distributions that characterize 
rather different environmental conditions. Although the 
obtained distribution is smooth, it is not very reliable. 
Therefore, a proper choice of the parameter d is impor- 
tant. For the three datasets considered in this article 
(i.e.. North American avifauna, NYSE, and transporta- 
tion network in Minnesota) we have chosen d = 1.5. We 
have checked that, in all the considered cases, the value is 
optimal. With d = 1.5, the shape of frequency distribu- 
tions is already sufficiently smooth and the obtained quo- 
tients, Pi{N)/ P2{N), do not change significantly com- 
pared with the smaller values of d (see Fig. [5]). 

For better understanding of data processing we have 
developed a web page where the data sets as well as a 
software for their analysis have been placed [isj . 



Appendix B: Supplementary equations 

1. Approximating function by its Poisson 
Transform 

Let us consider the Poisson transform defined as fol- 
lows: 



g{n) = P.T.[GiN),n]^J2GiN)- 



(Bl) 



N=0 



In this appendix, we show that the accuracy of the ap- 
proximation given by Eq. (jlOp . i.e. 



g{n) ~ G{n), 



(B2) 



may be evaluated on the basis the following theorem: 

Approximation theorem for the Poisson trans- 
form, as defined by Eq. (|B1|) . Let G{n) be the Pois- 
son transform ofG{N). Then we have 

j>i i=j+i 

where Q^^'^{N) corresponds to the i-th derivative ofQ{n) 
at n — N and 



k>0 



k 

k~J 



(B4) 



where 



k 

k-j 



stands for Stirling numbers of the first kind. 

One can show that the coefficients Cij satisfy the fol- 
lowing recurrence relation (ssj 



(i -I- l)ci+i j+i 



Ci-l, 



(B5) 
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FIG. 6: (Color online). Dependence of the quo- 
tients, Pi{N)/ P2{N), on the parameter d. a, Blue Jay 
{Cyanocitta cristata) populations in North America, (A'^i) = 
1.22 and {N2) = 3.71. b, Dynamics of the NYSE, (A^i) = 
6.62 X lO'^ and {N2) = 2.26 x 10^. c, Traffic intensity in 
Minnesota, (Ni) = 2.16 and {N2) = 8.59. 



with boundary conditions cqj — Sqj, where Sij stands 
for the Kronecker's delta, and cij — 0. Table lists 
some of these coefficients Cij. 

To prove the theorem we will use an analogous approx- 
imation theorem for the Poisson transform defined in the 
following manner 



J-„,(n) = P.T.* [F„,iN),n] - ^ ^,„(iV)- 



N 



N=0 



(B6) 

(Note that the Poisson transforms G{n) (|B1I) and J-,ri{n) 
(IB6p differ between each other. The former naturally 



emerges in many physical problems , e.g. in optics f3g 
and the science of complex networks [33, js^, whereas the 
latter is widely used in computer science and information 
theory.) The theorem states 'ss']: 

Approximation theorem for the Poisson trans- 
form, as defined by Eq. (|B6|) . Let J-m{n) be the 
Poisson transform of Fm{N). Then, for n = N/m, we 
have 



F™(7V)=J-,„(n)+5](l 



2j 



(B7) 

where J^m (n) corresponds to the i-th derivative ofJ-m(n) 
and the coefficients Ci,j are given by Eq. fi?^[ ). 

To prove the approximation theorem for the Poisson 
transform given by Eq. (jBll) we first put m = 1 in 
Eq. (|B7[) . Next, we note that the two transforms, G{n) 
(|Bip and J^i{n) l|B6|) represent the same function, i.e. 



when 



Fi{N) = G{N) 



m 



-N 



(B8) 



(B9) 



Then, replacing Fi{N) in Eq. dBTj with Eq. (HH, and us- 
ing Eq. (|B8p we immediately get Eq. (IB3p . This finishes 
the proof. 

Summarizing, approximating function G{n) by its 
Poisson transform Q (n) , cf. Eq. (|B2I) , is acceptable when 
the transform varies slowly enough. 



2. Analytical formula for the coefficients /„ 

In the main text of our paper, the parameters, 
/i,/2,---,/w stand for coefficients of the consecutive 
terms in the MacLaurin expansion ^ of the free en- 
ergy, F{fi) (ITSt . It means that the single coefficient, /„, 
corresponds to the 71th derivative of F{^) at = 



/„ = F(")(0). 



(BIO) 



Given Taylor's fluctuation scaling ([T]) and the 
fluctuation-dissipation relation ()13p . it is not diffi- 
cult to derive the general formula for /„. The flrst steps 
of this derivation are given below. 
First, one has 



/o = F{0). 



(Bll) 



The coefficient, /i, simply results from the fluctuation- 
dissipation relation 



/i = 



^Fi^l) 



= -(^>U=o. 



(B12) 
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TABLE II: Values of the coefficients c, 



Having the expression for /i,the coefficient, /2, can be 
calculated in the following manner 



/2 



.0 dfi J 



(B13) 



{{N')-{N)^)\,^o = {a{Nn,=o- 



The next coefficients /„ can be derived in a similar way. 
In particular, using Eq. (|B13|) one gets 



/3 = 



djajN)'') 



^=0 V dfi^ J 



(B14) 



Ai=0 



(-l)2a2fe(7V)(2''-i)U.o. 



Continuing these calculations, one can show that the gen- 
eral formula describing /„ for n > 2 can be written as 
follows 



d'^Fin) 



(B15) 



/i=0 



where 



n[("- 2)6 -(71 -3)], (B16) 



and (A^) is given by Eq. (IT4l) . 

Eqs. (|B15|) and (|B16|) have been used to prepare Fig.m 



3. Number of states characterizing the Poisson 
distribution 



The Poisson distribution. 



{N) 



N 



Nl 



(B17) 



fulfills Taylor's power law ([T]) with parameters a = b = 1. 
The mean value (I14p and the free energy (1151) correspond- 
ing with this distribution are given by 



(N) = Xe-^" 



(B18) 



and 



= {N)+Y = Xe-^ + y. 



(B19) 



where X and Y represent integration constants. In the 
following we calculate the number of states giN) char- 
acterizing the Poisson distribution in a twofold manner: 
first, taking advantage that we know the closed formula 
for the distribution (|B17|) . and second, using Eqs. 
and (fTTj) . In the second derivation one does not need to 
know the closed formula for the distribution (it is only 
known in a few cases, including the considered case of 
the Poisson distribution). The whole information about 
the system is taken from the free energy F{fi). 

Thus, putting Eq. (jBlSp into Eq. (|B17|) one gets the 
following form of the considered distribution 



iV! 



(B20) 



Similarly, inserting Eq. (jB19P into Eq. ([3]) one gets the 
equivalent formula for the distribution 



P{N;^l)=g{N)- 



-fj.N 



(B21) 



Comparing Eqs. (|B20)) and (|B2T|) one finds that Y — 
and g{N) corresponding to the Poisson distribution is 
given by 



x^ 
m ■ 



(B22) 



The expression (jB22l) for the number of states has a 
very simple interpretation. The numerator of g{N) cor- 
responds to the number of A^-variations (each of size N) 
from a set of size X. It means that the Poisson distri- 
bution arises in such systems, where each of N elements 
(birds, cars, etc.) can be found in one of X states, re- 
gardless of state of other elements. The denominator, 
A^!, of g{N) automatically assures us that elements of 
the considered system are indistinguishable. 

In the following, we show that the same formula for the 
number of states characterizing the Poisson distribution 
can be directly obtained form Eqs. (|8|- (|lip . To use the 
mentioned equations, one has to find the coefficients /„ 
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([T6)) in the MacLaurin expansion of the free energy 
(|B19[) . It is easy to see that the coefficients are given by 

/„ = i-irx. (B23) 

Inserting the coefficients into Eq. (jlip one gets 



71 I ^ ^ 



(B24) 



fc=i 



where S{n, k) — _B„^fc(l, 1, . . . , 1) stands for the Stirhng 
number of the second kind. To derive Eq. (|B24|) we have 
make use of the well-known properties of Bell polynomi- 
als [llj: the first one, 

n 

Bn{xi,X2j . . . ,a;„) = ^ Bn^k{xi,X2, ■ ■ .,Xn-k+l)j 
k=l 

(B25) 

stating that the nth complete Bell polynomial, i?„, is a 
sum of partial Bell polynomials, Bn^k, and the second 
one, 

Bn,k{abxi,ab'^X2, ■ ■ ■) = a''6"B„,fc(a;i, a;2, ■ • ■ )> (B26) 
resulting directly from definition of Bn,k^ 



Bn,k{xi,X2,- ■ ■) = 



ci!c2!...(l!)=i(2!)-2...-^i •••> 

(B27) 

where the summation takes place over all integers 
Ci,C2,C3, • • • > 0, such that ci -I- 2c2 -I- 3c3 + • • • = n 

and ci + C2 + C3 + ■ ■ ■ = k. 

The sum on the right hand side of Eq. (|B24p can be 
recognized as the right hand side of the so-called Dobihski 
formula lisl. 



1 t- 



N=0 



m 



(B28) 



fc=i 



that gives the nth moment of the Poisson distribution 
with expected value X. Now, putting Eq. (|B28p into 
Eq. (|B24| . and then comparing the resulting formula 
with Eq. ([5]) one finds the number of states g{N) given 
by the same formula, Eq. (IB22[) . as the one derived at 
the beginning of this subsection. 



4. Fitting of experimental frequency distributions 
in Fig. [2] 

Theoretical frequency distributions, P{N), describing 
systems with Taylors fiuctuation scaling have the follow- 
ing form (cf. Eqs. © and (^1]) ) 



PiN)=g{N)- 



(B29) 



To fit the experimental data with the formula, one has to 
find the correct value of the parameter, fj,, and also use 



the proper expressions for the free energy, F{ii), and the 
number of states, g{N). The value of /i can be calculated 
from the formula for the first moment of the distribu- 
tion (N) (mi). The function, F{fi), is given bv ([T5|). A 
problem, however, arises with g{N) underlying the con- 
sidered systems. The closed form expression for g{N) 
only exists in a few cases of the Taylor's parameter, b. 
Of course, in the limit of A'^ ^ 1, one can approximate 
g{N) with Eq. ([12]), but the numerical computation of 
the corresponding Bell polynomials, Bjq{fi, /2, . . . , /at), 
is extremely challenging and time-consuming. Therefore, 
in order to fit experimental frequency distributions in 
Fig. [2] with Eq. (|B29[) we have applied a numerical pro- 
cedure whose details are exposed below. 

Let us note that, although the formula for the fre- 
quency distribution, P{N) (jB29l) . does not have a closed- 
form expression for arbitrary values of > 1, the closed- 
form expression for generating function, P(s), of this dis- 
tribution does exist: 



P(s) 



JV=0 



P{N)s 



N 



,i^(M-lns)-_F(,x) 



(B30) 



Thus, the probability distribution, P{N), may be ob- 
tained from the series expansion of P(s) at s = 0, i.e.. 



PiN) 



1 

iV! 



dWp(.s) 



ds 



N 



(B31) 



s=0 



Unfortunately, here, due to the fact that In s and its con- 
secutive derivatives diverge at s = standard numer- 
ical procedures do not cope with Eq. (|B3ip . To over- 
come this problem, we have expanded the logarithm in 
(|B30|) into a power series up to the second order, i.e.. 

Ins ~ (s — 1) — -^^ — 2"^. Having a new generating func- 
tion. 



P(s) 



N=0 



P{N)s 



N 



(B32) 

we were able to calculate its consecutive derivatives at 
■3 = and also the corresponding approximated fre- 
quency distribution, P{N) (|B3ip . 

Fig. [2] in the main body of the article presents four 
experimental frequency distributions fitted with the ap- 
proximated probability distribution, P{N). In the Fig- 
ure, a direct fitting of the distribution describing stock 
market dynamics has been unsuccessful due to practical 
constraints related to impossibility of calculating deriva- 
tives of the desired order, 10^, for (IB31[) . Nevertheless, 
a theoretical frequency distribution characterizing NYSE 
has been obtained indirectly thanks to the scale- free char- 
acter of Taylor's power law. The procedure of the indirect 
fitting consisted of rescaling the parameter, N, with the 
factor, z, according to iV N/z. Taylor's law describing 
the rescaled systems has the following form 



(B33) 
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providing us with the new vahie of Taylor's parameter, 

a* = az''-^, (B34) 

which characterizes the frequency distribution in ()B3ip 
describing the rescaled parameter, N/z. 



poj . We also would like to point out that the idea of 
key factors is particularly interesting in relation to stock 
market dynamics and dynamics of other man-made sys- 
tems with Taylor's fluctuation scaling. Our derivations 
show that the systems possess a well-defined number of 
states and, therefore, a well-defined group of key factors. 



Appendix C: Supplementary discussion 

1. Key factor analysis 

The first remark relates to what we call key factor anal- 
ysis. Let us note that, although Eq. (fT2l) is the exact the- 
oretical result, which helps to understand the meaning of 
the number of states, in reality it may not be the sim- 
plest way to calculate the NoS. For example, in the case of 
three-dimensional, free electron g{E) describes the num- 
ber of states that are available to be occupied when the 
electron has energy E and is given by g{E) cx ^/E. The 
relation can be simply derived form the formula for the 
energy E{p^,py,p^) = {pi +pl +pl)/2m, which depends 
only on momentum components, Px, Py, and Pz- The con- 
stant energy surface corresponds to the sphere of radius, 

\JPx + Py + Pz VE. Therefore, the number of states 

with energy < E is given by T{E) oc E^^"^ , and the num- 
ber of states with energy E is g{E) = dr{E)/dE cx Ve. 

The above example shows that calculation of the num- 
ber of states, g{N), underlying Taylor's power law should 
start by specifying the key factors accounting for the 
quantity N in the considered systems (in the case of elec- 
tron, the key factors are momentum components). The 
idea behind our theoretical approach is that a set of such 
factors exists and significantly contributes to the varia- 
tion of N. An analogous idea underlies the key factor 
analysis (or the life table analysis) in population biology 
[sl]. Contrary to the classical approach to the problem 
proposed by biologists, however, our approach goes far 
beyond the analysis of the simplified key factor indices 



2. Dynamical processes leading to Taylor's law 

The second comment relates to dynamical processes 
leading to Taylor's law. It is important to understand 
that the proposed theoretical approach to the origins of 
the law provides static pictures of the phenomena. It 
does not explain details of the dynamical processes un- 
derlying the fluctuation scaling. The same is true for 
statistical mechanics, which reproduces static aspects of 
thermodynamics without dealing with dynamical pro- 
cesses such as the approach to equilibrium or dynamics 
in equilibrium. To investigate these points from the mi- 
croscopic dynamics is indeed a very difficult and poorly 
understood problem (see e.g. [HI). In most of cases, 
however, a rather satisfactory understanding of these is- 
sues is possible due to simple methods borrowed from the 
theory of stochastic processes and computational statisti- 
cal physics. For example, the detailed balance condition 
[4^ applied to the known microscopic models accounting 
for Taylor's law 23 25, 30] may help us to understand the 
meaning of the macroscopic parameter, /x cf. Eqs. (4) - 
(6). In this sense, our approach does not cancel the pre- 
vious models. It merely allows identification of models' 
important features that account for the law. Finally, let 
us note that the meaning of l//i is similar to the mean- 
ing of thermodynamic temperature, T. For physicists, 
it has taken a long time to understand the meaning of 
temperature and entropy. We believe that our approach 
will help to introduce and understand these concepts in 
other areas of science. 
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